Transition to zero resistance in a two dimensional electron gas driven with microwaves 
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High-mobility 2D electron systems in a perpendicular magnetic field exhibit zero resistance states (ZRS) 
when driven with microwave radiation. We study the nonequilibrium phase transition into this ZRS using phe- 
nomenological equations of motion to describe the current and density fluctuations. We focus on two models for 
the transition into a time-independent steady state. Model-I assumes rotational invariance, density conservation, 
and symmetry under shifting the density globally by a constant. This model is argued to describe physics on 
small length scales where the density does not vary appreciably from its mean. The ordered state that arises in 
this case breaks rotational invariance and consists of a uniform current and transverse Hall field. We discuss 
some properties of this state, such as stability to fluctuations and the appearance of a Goldstone mode associated 
with the continuous symmetry breaking. Using dynamical renormalization group techniques, we find that with 
short-range interactions this model can admit a continuous transition described by mean-field theory, whereas 
with long-range interactions the transition is driven first-order. Model-II, which assumes only rotational invari- 
ance and density conservation and is argued to be appropriate on longer length scales, is shown to predict a 
first-order transition with either short- or long-range interactions. We discuss implications for experiments, in- 
cluding scaling relations and a possible way to detect the Goldstone mode in the case of a continuous transition 
into the ZRS, as well as possible signatures of a first-order transition in larger samples. We also point out the 
connection of our work to the well-studied phenomenon of 'flocking' . 



PACS numbers: 73.40.-c,73.43.-f,78.67.-n 

I. INTRODUCTION 

High-mobility two-dimensional electron gases (2DEGs) 
subject to a perpendicular magnetic field exhibit novel physics 
when driven with microwave radiation. Zudov et al} first 
demonstrated that the longitudinal resistance develops dra- 
matic radiation-induced oscillations at low temperatures {T ^ 
\K) and low magnetic fields {B < IkG). These oscillations 
are periodic in 1/B, with the period set by the ratio of the 
microwave and cyclotron frequencies. In contrast, the Hall 
resistance is nearly unaffected by the microwaves', although 
small radiation-induced Hall oscillations have recently been 
observed^'^^. More spectacular is the observation, made in- 
dependently by Mani et alf" and Zudov et air', that in even 
higher-mobility samples the oscillations become sufficiently 
large that the minima of the resistance oscillations develop 
into zero resistance states — the measured resistance vanishes 
within experimental accuracy over a range of magnetic fields 
and radiation intensities. Subsequent experiments have con- 
firmed their results^'*-^ and also observed a similar effect in 
Corbino samples^, where zero-conductance states have been 
measured. 

On the theoretical front, several groups have caiTied out 
microscopic calculations of the resistance taking into account 
photon-assisted impurity scattering^iSiiiii^ii^iii and radiation- 
induced changes in the distribution functioni^. Although 
these mechanisms are quite different, both capture the resis- 
tance oscillations with the correct period and phase at low ra- 



diation intensity. At higher intensities, however, these calcula- 
tions predict a negative resistance in regions of magnetic field 
where the experiments find a zero resistance state. 

The missing ingredient needed to connect the microscopic 
theory with the experiments was the observation of Andreev 
et al. that a state characterized by a negative longitudinal re- 
sistance, quite independent of its microscopic origin, is un- 
stable to current fluctuations.'^ They argued that this instabil- 
ity leads to an inhomogeneous state where the system spon- 
taneously develops domains of cuiTent with magnitude jq, 
where jo coiTesponds to a vanishing longitudinal resistivity, 
i.e. pdUo) = 0. Applying an external current then merely 
reorganizes the domain sizes in order to accommodate the ad- 
ditional current, leading to zero measured resistance over a 
range of bias current as observed experimentally. In the ab- 
sence of an applied current, since the system has a large Hall 
resistance, this picture indicates that spontaneous cuiTent do- 
mains in the zero resistance state should reveal themselves 
through spontaneous Hall voltages transverse to the domains. 
Willett et alJ have indeed measured spontaneous voltages be- 
tween internal and external contacts with no applied current, 
which lends support to this idea. 

The experiments together with the microscopic calculations 
and phenomenological arguments provide strong evidence for 
the existence of a nonequilibrium phase transition from a nor- 
mal state with nonzero resistance to a zero resistance state 
whose detailed properties remain largely unexplored. In this 
paper, we attempt to gain an understanding of the nature of 
this transition, and to learn about the properties of the zero 
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resistance state. 



A. Strategy 

We begin with the observation that while the microscopic 
mechanism for how radiation induces the transition to a zero 
resistance state is a matter of some debate, this knowledge is 
not crucial for studying universal properties close to the phase 
transition. Indeed, in order to study the long-wavelength, low- 
frequency dynamics near the transition it is sufficient to iden- 
tify the appropriate hydrodynamic variables and construct the 
most general local equations of motion for them consistent 
with symmetries and conservation laws. The magnetic field, 
temperature, microwave radiation, and quantum effects will 
determine the various parameters of this theory; these may 
be calculated in principle from a microscopic approach, but 
we do not attempt to do this here. Our idea will be to view the 
equations of motion as a non-equilibrium analogue of Landau- 
Ginzburg-Wilson theory. We will use them to study universal 
physics near the phase transition, going beyond mean field 
theory by including nonlinearities and fluctuations within a 
renormalization group framework. 

In the vicinity of the transition into the zero resistance state 
in the 2DEG, the relevant hydrodynamic degrees of freedom 
are the current density j(r, t) and the charge density n(r, t), 
which are constrained by a continuity equation, 

dn 

^ + V.j^O, (1) 

that enforces local charge conservation. 

The dynamics of the current density j (r , t) is governed by a 
nonequilibrium equation of motion (akin to the Navier-Stokes 
equation) for a 2D charged fluid in a perpendicular magnetic 
field. Because of the nonequilibrium nature of the system 
(microwave-driven 2D electron liquid) the equation for j in- 
cludes non-conservative forces, i.e., those not derivable from a 
free energy functional. Hence, the generic symmetry-allowed 
form of the equation for the current density is only restricted 
by the translational and rotational invariances in the plane. 
Keeping the leading order (at long length and time scales) 
terms in powers of the charge and current densities and their 
gradients leads to 

LOo'dfi + 9ij = -rj-u|j|2j + ,yiV2j + ?72V(V-j) 

- Mi ■ v)j - i^2Vj2 - i/3(v • j)j + 7i$v$ 

+ 72«'j + 73^-2 X j + C + • ■ • . (2) 

As we will discuss in more detail below, terms appearing on 
the right-hand side of the above equation are forces that de- 
termine the local acceleration (9tj(r, t)) of the electron fluid, 
each having a simple physical interpretation. The r and u 
terms are the linear and nonlinear longitudinal resistivities 
(frictional drag forces on the electron fluid). The rji terms 
describe viscous forces associated with a nonuniform flow 
and the ujc term is the Lorentz force on the charged moving 



electron fluid. The i/^ terms are convective-like nonlineari- 
ties, where the absence of Galilean invariance permits more 
general types of convective terms 1^2 and 1^3 in addition to the 
conventional i/i term, with generic (symmetry unrestricted) 
values of these couplings. 

Here, the potential $ = is determined by the density 
via 

$ = ^ y(r-r')n(r'). (3) 

For long-range interactions, V{r — r') ^ l/|r — r'| is the 
Coulomb potential. For a screened interaction, we can set 
V{r - r') « S(r - r'), so that $ « 71. With this, the fi 
and 7i terms incorporate Fick's law (diffusion down a local 
chemical potential gradient), with the latter accounting for a 
density-dependent diffusion coefficient. Similarly, 72 and 73 
account for the lowest order density-dependence of the linear 
resistivity and the Lorentz force. 

In addition, we have included in Eq. Q a zero-mean white 
noise force ( with a correlator 

(r (r, t)C^{r\t')) = 2gS'^^'6{r - r')S{t - t'). (4) 

Apart from thermal noise, this incorporates the effect of mi- 
croscopic fluctuations that arise from the coarse graining im- 
plicit in our formulation. Since we are dealing with a system 
far from equilibrium, the strength g of the noise is not fixed 
by the fluctuation-dissipation relation, but is an independent 
quantity. 

Focusing on the terms dt] = — rj + in Eq. (|2ji, it is clear 
that (i) for large positive values of r, the zero current state 
is stable and current fluctuations decay exponentially, while 
(ii) for large negative values of r, current fluctuations grow 
exponentially and the zero current state is unstable. Thus, as r 
changes from positive to negative (in the experiments tunable 
by a microwave power and/or frequency), Eqs. ( Ill4t describe 
the phase transition from a conventional resistive state for r > 
to a nonequilibrium steady state with spontaneous currents 
for r < 0. 

As we will argue in Section II, this set of equations can po- 
tentially describe various types of current and density order- 
ing, including circulating current states and domain patterns 
of current and density. In this paper our main focus will be 
on the nature of the transition into a time-independent steady 
state with possible density and current domains since, given 
the observations of Willett et ai, this appears to be relevant 
to the 2DEG experiments. We defer to future work questions 
regarding the detailed nature of the ordered state in this case, 
as well as a study of the phase transition into the circulating 
state. 

A quite different theoretical motivation for studying this 
problem arises from the observation that the current and den- 
sity evolution equations studied here reduce, for B = and 
short range interactions, to the continuum equations used to 
investigate the problem of 'flocking' . "^ '"^ In that case, the sys- 
tem has been shown to develop an expectation value for the 
particle current, thus spontaneously breaking the continuous 
rotational symmetry even in two spatial dimensions. This is 
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particularly striking since the Mermin-Wagner theorem ■ for- 
bids such symmetry breaking in d = 2 for classical equilib- 
rium systems. This 'violation' was identified as arising from 
nonlinear convective terms which are only allowed in non- 
equilibrium systems, and turn out to be relevant for this prob- 
lem in dimensions d < 4. Much is known about the univer- 
sal dynamics in the flocking state in d = 2, but the nature 
of the phase transition into this state has not been addressed 
analytically. The question we study is equivalent to asking: 
What is the fate of the flocking transition and the flocking 
state in two dimensions in the presence of a magnetic field 
that breaks time-reversal symmetry? As we show, one can 
make more progress in this modified problem. This 'flocking' 
point of view is also useful for carrying out numerical simu- 
lations, since many simple particle models for flocking have 
been studied in the absence of a magnetic field and can be 
adapted to our problem, although we do not pursue this here. 

B. Summary of the paper 

We begin in Section II by showing how some terms in 
Eq. (13 can be related to the full nonlinear resistivity. We do 
this by formally expanding the relation 

E(k,cj) X £ (5) 

at low frequency and wavevector, and for small current and 
potential fluctuations. Here pD and pn represent the diag- 
onal and Hall resistivities, and the electric field E is deter- 
mined via the electrostatic potential, i.e. E — — V$. Upon 
Fourier transforming back to real space, one can arrive at an 
equation with a form similar to Eq. This proves to be 
a useful exercise since we can then relate different possible 
forms of the frequency and wavevector dependent resistivity 
in the presence of microwaves to the model parameters ap- 
pearing in Eq. (|2j and therefore to the kinds of ordered states 
which might emerge from our description. More importantly, 
this helps us to identify the correct set of critical modes near 
the phase transition into these putative ordered states. Specif- 
ically, we show that if the resistivity is an increasing function 
of frequency at low frequency, so that the zero resistance state 
is achieved when the DC resistance first goes negative, then 
a time-independent steady state with inhomogeneous density 
would result. The only critical mode near the transition into 
this state involves density fluctuations accompanied by current 
fluctuations that balance the Lorentz force. Since the current 
and density are tied to one another in this mode, one can re- 
express current fluctuations in terms of the density. Inserting 
the resulting expression into the continuity equation results in 
an equation of motion involving only the density at the critical 
point. 

This equation of motion for the density at the critical point 
depends on terms involving the absolute magnitude of the 
density, as well as terms that depend only on density gradi- 
ents. We warm up in Section III by analyzing a model which 
neglects terms that depend on the absolute magnitude of the 
density, and is instead invariant under shifting the density by 
a constant. This model is expected to describe physics on 



short length scales where the density does not vary apprecia- 
bly from its mean so that such terms can be safely ignored. 
In the 'ordered phase' of this model the system develops a 
uniform current with a transverse density gradient that bal- 
ances the Lorentz force. We show that this state is stable to 
small fluctuations, and discuss the 'Goldstone mode' associ- 
ated with the spontaneously broken rotational symmetry. The 
ordered state described by this model is argued to be relevant 
for the experiments at short length scales L < Ld, where Ld 
is estimated to be roughly 1mm, comparable to sample sizes 
used in the experiments. We then turn to the critical properties 
of this model, considering both short- and long-range interac- 
tions. With short-range interactions, we show that the upper 
critical dimension is due = 2. In this case, we use dynami- 
cal renormalization group calculations to demonstrate that the 
Gaussian fixed point has a finite-volume basin of attraction; 
hence, a finite fraction of initial nonlinear couplings all flow 
to zero upon renormalization. In such cases, the transition is 
continuous and governed to a good approximation by mean 
field theory. Various scaling relations should hold near the 
transition in this regime. For instance, at fixed magnetic field 
strength and in the absence of an applied voltage, below the 
transition, the spontaneous current jo should scale with the 
microwave power P as 

jo{P)-\P-Pcf, (6) 

where (3 > and Pc is the critical microwave power at which 
the longitudinal resistance first vanishes. Approaching the 
transition from the resistive side, with P < P^, we expect a 
universal scaling relation to hold between the imposed current 
i and the induced longitudinal electric field, 

E) ^ (P, - Pff (E^/'{Pe - P)-P) , (7) 

where f{x) is a scaling function with the properties that 
f{x) as a; ^ 0, and f{x) ~ a; as x — > cx). The be- 

havior as a: ^ recovers the linear response result, j oc E, in 
the resistive phase, with a resistivity ^ [Pc — P)f'^^^^\ The 
behavior for x — > oo leads to a universal longitudinal nonlin- 
ear IV characteristic 

jiPcE) E'^' (8) 

at the transition, P = Pc, with mean-field value of (5 = 3 for 
a current-biased experiment. We then show that long-range 
interactions appear to drive the transition first-order. The ex- 
perimental signatures of the Goldstone mode in the 'ordered 
phase' and the mean field transition with short range interac- 
tions are qualitatively discussed in Section V. 

In Section IV we analyze the phase transition in the more 
general model, where terms that depend on the absolute mag- 
nitude of the density are taken into account. These terms, 
which become important on length scales L > Lc2, are ar- 
gued to drive the transition first-order with either short- or 
long-range interactions based on renormalization group cal- 
culations. We derive an expression for Lc2 that depends on 
the density- and wavevector-dependent resistivity, and sug- 
gest that microscopic calculations may be used to estimate this 
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length. Experimental consequences of the first-order phase 
transition for sample sizes larger than ic2 are briefly noted in 
Section V. 



II. DERIVING AND SIMPLIFYING THE EQUATIONS OF 
MOTION 




A. "Microscopic derivation" of equations of motion 



Before we turn to the analysis of the phases and transitions 
described by the set of Eqs. ill4\ . let us consider a derivation 
of some terms in the equation of motion for j in Eq. (|2j- 

We begin with the linear response relation 

E"(k, iv) = p%^ik, u;)fik, Lu) + pH{k, c.)e"^j^(k, c.) (9) 

where and pn represent the diagonal and Hall resistivi- 
ties, the electric field E is determined from the electrostatic 
potential via E = — V<i>, and e"^ is the antisymmetric tensor. 

We know from experiments that ph{0, 0) « B/ (ne) even 
in the presence of microwave radiation. We are interested in 
the case where the dissipative part of the microscopic diago- 
nal resistivity becomes negative. With increasing microwave 
intensity, this would first happen at some particular wavevec- 
tor and frequency (K, il). Two specific cases for the behavior 
of pdO^, are illustrated in Fig. [J 

If K, 57 are small, we can access the resistivity mini- 
mum shown in Fig.[2by expanding /9£)(k, w) = pi(k, w) + 
ip2{k, uj) in a Taylor series as: 



a/3 

Pd 



ga,3 
772k"k'3 



(10) 



where the frequency derivatives and coefficients 771.2 are eval- 
uated at (k = 0,cj = 0). Using this expansion inside Eq. 
(|9j, and assuming that the Hall resistivity is independent of 
wavevector and frequency in the regime of interest, we find 
the following relations between the coefficients in Eq. Q and 
the microscopic linear response resistivity. 



r 


= Pi(0,0)/G 


(11) 




1 d^Pi 
2G duj^ 


(12) 




= l/G 


(13) 




= m/G 


(14) 


m 


= mIG 


(15) 




= Ph(0,0)/G, 


(16) 



where G = —(dp2/duj). 

We can similarly match some of the non-linear terms in 
Eq. (O as follows. Let us take k = 0, = and consider 
the nonlinear resistivity that depends in general on the local 
potential and the current magnitude, namely. 



PHif,^) 



= '5"^[pi(0,0) 
= Ph{0,0) + UHi 



2 ^-3$ + . 



(17) 
(18) 



FIG. 1: Schematic behavior of the real part of po when (a) ljo < 
and (b) uo > 0- 



Using this expansion and comparing with the non-linear terms 
in Eq. we find 



u — u/G 

72 = 72/G 

73 = 73/G. 



(19) 
(20) 
(21) 



As we shall see below, the type of ordering expected to emerge 
from our description depends on the frequency and wavevec- 
tor dependence of the resistivity - measuring these in the dis- 
ordered phase close to the transition would offer clues to the 
nature of the zero resistance state. 



Identifying critical modes and simplifying the equations of 
motion 



On general grounds, one would expect that the type of or- 
der that develops near the transition should depend on where 
the minimum of pi occurs in (/c, lo) space. The resistance will 
in general depend on both k and uj, and close to the transi- 
tion will only be negative in a small region of frequencies and 
wavevectors about the minimum. Modes away from the mini- 
mum remain stable. Two possibilities for where this minimum 
occurs as a function of frequency are sketched in Fig.^ If the 
minimum occurs at zero frequency as in Fig.^a), then zero- 
frequency modes that become critical at the transition should 
give rise to a time-independent ordered state (e.g., static do- 
mains of current). If on the other hand the minimum occurs 
at a non-zero frequency as in Fig.^b), then finite-frequency 
modes should give rise to a state ordered at finite frequency 
(e.g., circulating currents). In either case, the wavevector at 
which the minimum occurs would determine the wavevector 
at which the system orders. Thus, the signs of 771 and c^o, 
which determine whether the minimum of pi occurs at zero 
(or nonzero) wavevector and frequency, should play an im- 
portant role in the ordering. 

To make this more concrete, let us consider the mode struc- 
ture in the disordered state, where j and n represent fluctua- 
tions about a stable zero-current state. We start with the case 
r/i > and focus on wavevectors k ^ Q since the resistiv- 
ity is minimized when fc = 0. The modes obtained from the 
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linearized equations of motion are given by 



/ 2r 

±L0Al + — 
UJ[) J \ Wo 



(22) 
(23) 



where V{k) is the Fourier-transform of the interaction poten- 
tial V{r). Equation (I22> is only written out to order oj'^^ for 
simplicity. We only want to consider here the effect of adding 
a small frequency dependence to the resistivity, so the exact 
expression is not important. The modes in Eq. \22\ correspond 
to current fluctuations that circulate due to the magnetic field 
as they dissipate. The associated density fluctuations for these 
modes vanish in the fc — > limit. Equation ( I23> represents a 
diffusive mode involving both current and density fluctuations 
that survive in the /c — > limit. These current fluctuations are 
undeflected by the magnetic field because the Lorentz force is 
balanced by an electric field set up by density fluctuations. 

In order for the zero-current state to be stable, the imaginary 
part of these frequencies must be negative so that fluctuations 
are damped exponentially in time. For the diffusive mode, 
stability requires r > 0. The circulating current modes are 
stable when r > loI/ujq, assuming cjg 3> for simplicity. 
Violation of either inequality renders the zero-current state of 
the system unstable to current fluctuations. Since r ^ pi, this 
instability occurs approximately where the longitudinal resis- 
tivity changes sign, consistent with the findings of Andreev et 

As the longitudinal resistance tends to zero and the ordered 
state is approached, the circulating current modes become 
critical before the diffusive mode if ujq > 0. (Note that since 
these modes propagate at a finite frequency, this is consistent 
with the above discussion on the frequency-dependence of the 
longitudinal resistivity.) Once the cumulating current modes 
become unstable, the system should undergo a transition into 
an ordered state where circulating currents spontaneously de- 
velop but the density remains uniform. If ojq < 0, however, 
the diffusive mode becomes critical while the circulating cur- 
rent modes remain damped. In this case one would expect the 
system to undergo a transition into a phase with nonuniform 
density and spontaneous currents ordered at zero wavevector 
A distinguishing characteristic of the latter phase would be 
the development of voltages resulting from the nonuniform 
density. Since spontaneous voltages in the absence of a net 
current have indeed been observed in the ordered state, the 
case ojQ < seems to be the experimentally relevant one. We 
consequently focus on the transition into the density-ordered 
state and leave an analysis of the circulating-current state to 
future studies. 

These same ideas can be applied to the case rji < 0, where 
the resistivity is minimized at finite wavevector. Assuming 
773 > 0, one is then interested in wavevectors with magnitude 
close to fco = (1^71 l/S^ys)^^^, corresponding to the resistivity 
minimum. Since we can no longer perturb in k, we cannot in 
general write down simple expressions for the modes in the 
disordered state. We will therefore focus on the point where 



the resistance at zero frequency and k = ko drops to zero 
since this simplifies the mode structure. (This happens when 
r = 77^/4773.) A critical diffusive mode then emerges whose 
frequency is given to lowest order by 



UJD 



-iliklV{kn 
^fcgF(fco) + w 



(24) 



where 5k = |k| — fep. The circulating current modes to lowest 
order are 



UJ± 



-iklfi/2 ± ^^k^f,y4 + ^ilViko)+ul 



(25) 



where ?7 = 772 + 774^0- We have set cOq^ — here since the 
modes already do not become critical simultaneously. In the 
limit where the ujc term is dominant in Eq. ( I25t . the square 
root is positive. We will assume that 77 > so that these 
modes remain damped when the diffusive mode becomes crit- 
ical since this appears to be the experimentally relevant situ- 
ation. As the resistance decreases further, one would expect 
the diffusive mode to give rise to a time-independent state with 
nonuniform density ordered at wavevector fco- 

It follows from the preceding discussion that only the dif- 
fusive mode should be important for describing the transition 
into a nonuniform density phase ordered at either zero or finite 
wavevector. Since the circulating current modes have a finite 
damping rate when the diffusive mode becomes critical, they 
can be neglected provided we focus on frequencies smaller 
than their decay rate. This provides a large simplification in 
that it allows us to eliminate the currents altogether and ob- 
tain a theory in terms of the density alone. Physically, this is 
possible because at long times scales the current and density 
fluctuations are dominated by a diffusive mode characterized 
by a gradient of the density fluctuations that just balances the 
Lorentz force associated with the current fluctuations. One 
would thus expect to be able to write the 'fast' current in terms 
of the 'slow' density. This can be done by dropping the time 
derivatives on the left-hand side of Eq. (|2ji compared to ujc and 
then solving order by order for the current as a function of the 
density. Inserting the resulting expression into the continuity 
equation yields a decoupled equation of motion for the den- 
sity alone. The transition within this simplified description 
of the system will be analyzed in Sections III and IV for the 
case of zero wavevector ordering; finite-wavevector ordering 
is briefly mentioned in Section V but will not be studied in 
detail here. 



III. TRANSITION TO DENSITY-ORDERED STATE AT 
ZERO WAVEVECTOR WITH $ ^ <1> + const SYMMETRY 

When 771 > so that the resistivity is minimized at zero 
wavevector, we saw in the previous section that the zero- 
current state of the system becomes unstable when r < 0. 
Identifying the precise ordered state that develops in this 
regime is complicated by the presence of nonlinear terms 
in Eq. (|2ji involving the magnitude of the potential <E>. If 
one ignores such terms by manually imposing the symmetry 
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$ — * $ + const, then a simple ordered state emerges, namely 
a state with a uniform current and a transverse electric field 
that balances the Lorentz force. We will begin this section 
by discussing some mean-field properties of this ordered state 
and then analyze the transition to this state using dynamical 
renormalization group techniques. Our motivation for study- 
ing this simplified model is as follows. First, it is the simplest 
model that one can construct that captures the instability that 
occurs when the resistance becomes negative. Second, we ex- 
pect this model to be appropriate for describing physics on 
length scales where terms involving the magnitude of $ play 
a relatively unimportant role. This will be quantified below. 
Third, understanding the properties of the transition in this 
minimal model will allow us to better understand the effects 
of adding in terms that violate the <i> ^ $ + const symmetry, 
which will be done in Section IV. 



A. Ordered state and linearized tlieor y of fluctuations 

When r < 0, the ordered state within a model with $ 
$ + const symmetry consists of a uniform current 



imaginary parts are given by 



jo = vHA 



(26) 



where the direction x is spontaneously picked out. To balance 
the associated Lorentz force requires an electric field given by 



Eq 



-V$o = (wc/M)jo X 



(27) 



We have assumed here that jo is small in some sense so that, 
for instance, terms in the equation of motion proportional to 
can be neglected compared to the term. To further 
simplify things, terms such as (V<I')^j that would arise from 
expanding the longitudinal resistivity to higher order in the 
potential have also been neglected. Their presence only alters 
quantitative properties of the ordered state. For instance, a 
uniform current still develops, but with a modified magnitude. 

To establish a connection with the experiments, note that 
the longitudinal resistance at zero wavevector and frequency 



is proportional to 



- u\y\ (neglecting higher-order terms 



in j and V$). The spontaneous current jo therefore corre- 
sponds to a vanishing longitudinal resistance as seen experi- 
mentally. This is also consistent with the results of Andreev 
et al}'' that show that a stable state must have spontaneous 
currents corresponding to a vanishing longitudinal resistivity. 

To analyze the stability of the ordered state, we consider 
fluctuations about the uniform current state by writing j = 
jo + (5j and n — n^ + 6n, where jo is given in Eq. ( I26l l and 
no corresponds to the potential $o given in Eq. J27> . In the 
linearized equations of motion for Sj and Sn, there are two 
damped modes in the fc ^ limit with frequencies 



(28) 



Since the ordered state breaks rotational symmetry, there is 
also a Goldstone mode with frequency lliq whose real and 



ImujG 



(29) 

r;ifc|)t/(fc), (30) 



where k± and fc|| are the components of k perpendicular and 
parallel to jo, respectively. Note that the damping within this 
mode is anisotropic. In particular, fluctuations with wavevec- 
tor parallel to jo, which produce long-wavelength variations in 
the direction of jo, relax much more slowly than fluctuations 
with wavevector perpendicular to jo. 

To see if the ordered state is stable to fluctuations, one 
needs to calculate the mean-squared fluctuations of dj{r,t) 
and 6n{r, t), averaged over the noise. A divergence of either 
of these quantities would signal the destruction of the ordered 
state. We compute these quantities within the linearized the- 
ory, focusing only on fluctuations arising from the Goldstone 
mode for simplicity. (Long-wavelength fluctuations arising 
from the ut'^ modes will be finite since they have a nonzero 
damping rate as fc ^ 0; consequently, these modes can be 
neglected.) Denoting the current fluctuations parallel and per- 
pendicular to jo by (5j|| and (5j^, respectively, we find 
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liV{k)k'i{Coie - 



ak\ + T^ifc* 



g f i2V{k)[u^Pkl - a^{Col + a^)ki 
ak\^ + rjik^ 



./k 



(31) 



(32) 



where a = 2itjo and g is the noise strength. Equation 
(13 1> is obviously finite with either short-range interactions 
{V{k) const) or long-range interactions {V{k) ^ 1/fc) 
since the integrand itself is not infrared divergent. With long- 
range interactions, the integrand in Eq. ( I32t is infrared diver- 
gent. However, this divergence is integrable in 2D, leading to 
finite transverse current fluctuations. The mean-squared den- 
sity fluctuations are given by 



{Sn^iv,t)) 
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■ 2aujrk\\k I 



nV{k){akl+riik^) 



(33) 



Again, since the infrared divergence in Eq. ( I33> is integrable 
with either short- or long-range interactions, the density fluc- 
tuations are also finite. Hence we conclude that, for suffi- 
ciently low noise, g, the spontaneous current-carrying state is 
stable to current and density fluctuations with either short- or 
long-range interactions. 

The state characterized by Eqs. i26\ and ( I27l i can clearly not 
exist in arbitrarily large samples since the density would even- 
tually become negative on one side of the sample. For a given 
spontaneous current jo, one can estimate the maximum sam- 
ple length Lci below which this is a sensible ordered state by 
finding how large the sample can be before the density change 
becomes comparable to the mean density. We do this by as- 
suming that the electron-electron interactions are screened so 
that the electric field is given by the gradient in the electro- 
chemical potential /i. If A/i is the change in /i between the 
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edges of a sample of length L, then the magnitude of the elec- 
tric field is E — A/i/eL, where e is the electron charge. Re- 
gions of density variation comparable to the mean density will 
appear if A/i ~ Ep, where Ep is the Fermi energy. Setting 
A^ = Ep and using Eq. ( I27> . we get 

Lci - (34) 
epHjo 

where we have identified cDc//^ = Ph- 

Note that as the mean-field critical point is approached, 
jo — * and so Ld diverges. One might therefore be tempted 
to conclude that the model with $ ^ $ + const symmetry 
correctly describes the physics at the transition at all length 
scales. We stress that this is not necessarily the case. In com- 
puting Lci, we have only demanded that no unphysical fea- 
tures such as negative density arise in this minimal model. 
What we have not done is compute the characteristic length 
(which can be smaller than Ld above) below which terms 
that depend on the magnitude of $ play a negligible role. We 
will elaborate further on this in the following subsection. 

We now estimate Ld using parameters measured by Willett 
et a/fiin order to get a feel for this length scale. In their exper- 
iments, carried out on GaAs/AlGaAs samples, the density is 
n « 2 X lO^^cm^^, from which we estimate Ep ^ 5meV. In 
a 20GHz microwave field the primary zero resistance state oc- 
curs at i? w 0.4kG, where pn = B/ne « 125il. From spon- 
taneous voltages that develop in this zero resistance region, 
they estimate a spontaneous current of roughly 5/iA flow- 
ing between the center and edge in square samples of length 
0.4mm. Assuming a single domain between these contacts, 
we find jo ~ 25/iAmm^^. Putting these parameters together, 
we estimate the critical length for to be Ld ~ 1mm. In sam- 
ples with dimension larger than Ld, terms in the equation of 
motion involving the magnitude of $ must be taken into ac- 
count to produce a sensible ordered state. Such terms would 
prevent the density from becoming arbitrarily large and nega- 
tive, and would lead to inhomogeneous currents and densities. 
Determining the corresponding current-carrying ordered state 
on these longer length scales is an interesting problem that we 
do not address here. 



B. Transition with short-range interactions 

Having discussed an example of a stable ordered state that 
arises from a model with <i> ^ $ + const symmetry when 
r < 0, we now turn to the critical properties of system at the 
phase transition. We begin with the simplest case of short- 
range interactions. 

As discussed in Section II, our analysis is greatly simpli- 
fied by assuming that near the critical point circulating cur- 
rent modes remain damped while the diffusive mode become 
critical. Focusing only on the diffusive mode enables us to 
write the current in terms of the density and use the continuity 
equation to write an effective theory involving only the den- 
sity. Imposing the symmetry $—>(!> + const, the resulting 



equation of motion takes the form 

^ dtn- fV2$ + DV*$ - V • C' + AiV • (V^^V*) 
+ Asz • (V$ X V^$) + A3V • ((z X V$) • V(z X V$)) 
+ /iiV • (V$(V$)2) + Ai2z • (V$ X V(V$)2), (35) 

where D > and C is a Gaussian noise source with variance 
2g'. The transition occurs at f = (within mean-field theory), 
so we will take f — from now on. Since we are considering 
short-range interactions here, $[71] = n. Quartic and higher 
order terms in V<I> have been neglected since they contain at 
least five gradient operators and are therefore irrelevant in two 
dimensions. 

The upper critical dimension for this model is due — 2 for 
the case of short-range interactions. We use dynamical renor- 
malization group techniques^'-^^ to deduce whether the non- 
linearities appearing in Eq. i35\ are marginally relevant or ir- 
relevant in two dimensions. This procedure is facilitated by 
the use of the Martin-Siggia-Rose (MSR) formalism.^^. The 
essence of this formalism is that one introduces a 'partition 
function' Z that is useful for obtaining various correlation 
functions, namely 

Z = y VnS{dtn - fV^* + DV^<^> -V ■(+■■■)■ (36) 

This imposes the equation of motion as a constraint on all pos- 
sible spacetime 'trajectories' of n(r, t). The ellipsis indicates 
all nonlinearities appearing in Eq. i35\ . This functional delta- 
function constraint is implemented through an auxiliary field 
n(r, t) so that Z can be written as 

Z = j PnPne*^["'"l (37) 

S ^ [ n[at7i-fV2$ + DV*$- V-C' + ---],(38) 

Jr,t 

where constants have been absorbed into the integration mea- 
sure for n. A useful feature of this method is that the noise 
averaging can now be easily performed, with the result that 

f n[-V • C'] ^ ig' f (Vn)2 (39) 

Jr,t Jr.t 

in the 'action' S. This leads to an MSR action expressed in 
terms of the fields n, h and no noise terms. One can then im- 
plement the renormalization group transformation using stan- 
dard field theory techniques as follows. First, the action S is 
written in Fourier space with an ultraviolet cutoff A reflect- 
ing the coarse-graining of the fields. One then integrates out 
fields with wavevectors q such that A/s < g < A, where 
s > 1. This results in an effective action with a reduced cutoff 
A/s. To restore the initial cutoff, one then rescales wavevec- 



tors, frequencies, and the fields according to 

k' = sk (40) 

LU' = S^UJ (41) 

n'(k',w') = s"Xn(k,w) (42) 

n'(k',w') = s^*n(k,w). (43) 
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By setting s = 1 + d£, one then obtains differential recursion 
relations that specify how the effective coupling constants for 
the long-scale degrees of freedom "flow" as short-scale de- 
grees of freedom are integrated out. In the present paper these 
recursion relations will be calculated to one-loop order 

In anticipation of finding a stable Gaussian fixed point, we 
choose the rescaling exponents to take on their mean-field val- 
ues: z — A, X = and x = 4. (Since there is no small param- 
eter at our disposal, the only possible controlled fixed point 
must be Gaussian.) These exponents keep the noise strength 
g' fixed under renormalization since diagrammatic corrections 
to g' vanish at one-loop order To simplify the flow equations 
for the remaining coupling constants, we define the following 
dimensionless parameters: 



Ai = (a/i?)A3(4Ai -3A3) 

A2 = -(a/i:')A3A2 

A3 = {a/D)Xl 

A4 — — a/ii 

A5 = a^2, 



(44) 



where a ~ g/AnD^. The flow equations in terms of these 
parameters are 





-AiD 
2 






(45) 


diAi = 


-(-Ai 


+ 13A4)Ai - 3A3A4 


+ 4A2A5 


(46) 


diA2 = 


-(-Ai 


+ 7A4)A2 + i(3A3- 


- 7Ai)A5 


(47) 


diAs 




+ 12A4)A3 




(48) 


diAi = 


-(Ai + 


9A4)A4 + A^ 




(49) 


deA5 = 


-(Ai + 


10A4)A5. 




(50) 



At this point we would like to identify the basin of attrac- 
tion for the Gaussian fixed point under consideration. That 
is, for a given set of initial conditions for Aj, we would 
like to know whether these parameters all flow to zero as 
i ^ 00. While it is straightforward to check this numeri- 
cally, it is difficult to draw general conclusions either analyt- 
ically or from the numerics due to the five-dimensional pa- 
rameter space and the fact that Eqs. ( I46> through i50\ are all 
coupled. In the subspace with A3 = 0, one can show analyti- 
cally that the Gaussian fixed point is stable to all perturbations 
within that subspace. In the full parameter space with A3 7^ 0, 
we have shown that a finite-volume region of initial condi- 
tions corresponds to stable trajectories where each A; flows to 
zero. The asymptotic solution for such trajectories is given by 
Ai - (2/11)^-1, A2 - (5c2/242)£-io/", A3 - ci^-i^/ii^ 
A4 - (1/11)^-1, and A5 - (l/c2)^~^2/", where ci,2 are 
arbitrary constants. One can verify the stability of these flows 
by perturbing around this solution. According to Eq. ( I45> . the 
subdiffusion constant grows asymptotically as _D ~ Doi^^^^ 
along these trajectories, where Dq is a constant. The asymp- 
totic behavior of the original coupling constants is given by 
Ai ^ ^-2/11, A2 - £-1/", A3 ~ £-6/11, ^1 - £-9/", 
and fi2 ~ demonstrating marginal irrelevance of all 



nonlinearities, and therefore the stability of the Gaussian fixed 
point. 

Flows that terminate along the above asymptotic trajecto- 
ries correspond to marginally irrelevant couplings that reside 
in the basin of attraction for the Gaussian fixed point. In 
such cases mean-field theory should be a reasonable starting 
point for analyzing the transition to the density-ordered phase 
within this model. In particular, as we saw in the previous 
subsection mean-field theory predicts a continuous transition 
since the spontaneous current and the associated density gra- 
dient develop smoothly from zero as r becomes negative (see 
Eqs. ( I26t and ( I27» . Another important mean-field prediction 
we can make is that at the transition there is a single subdiffu- 
sive mode for density fluctuations with frequency 



-iDk^ 



(51) 



This slow relaxation of long-wavelength fluctuations should 
be accompanied by large voltage fluctuations near the transi- 
tion. Equal-time density-density correlations, which should 
mimic voltage correlations, are given within mean-field the- 
ory by 



(n(r,t)n(0,i)) 



9_ 
D 



fc2 



(52) 



The integral diverges logarithmically at small k. To regulate 
the integral, we restrict the range of integration to 27r/L < 
k < A, where L is the system size. In the limit rA 3> 1 and 
r/L <C 1, we obtain 



(n(r,i)n(0,t)} 



9' 



2ttD 



ln(i/r). 



(53) 



Equal-time current-current correlations in mean-field theory 
are given by 



(j(r,i)j(0,i)) 



g'^iA^ Ji(rA) 
27rD{:;^ rA ' 



(54) 



where Ji (x) is a Bessel function of the first kind. Note that the 
current becomes (5-function correlated in the limit that A 
00. Since the interactions are only marginally UTelevant, they 
will give rise to logarithmic corrections to these correlation 
functions, which will not be computed here. 

So far we have focused only on the case where the cou- 
pling constants flow to zero upon renormalization. Even out- 
side of this marginally stable region, the couplings are still 
only marginally relevant, and therefore grow only logarithmi- 
cally with length scale. In fact, we have seen numerically that 
many trajectories that initially flow toward the Gaussian fixed 
point eventually diverge from it, but do so only after many 
renormalization group iterations. In these instances, it may 
be very difficult to resolve deviations from mean-field theory 
either numerically or experimentally and the transition may 
appear continuous even in the presence of the marginally rel- 
evant couplings. 
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C. Is this model valid near the transition? 



D. Transition with long-range interactions 



We will now discuss when the model with <!> ^ $ + const 
symmetry and short-range interactions is expected to be ap- 
propriate for describing the physics at the transition. Consider 
adding the term AV • ($V<1>) to Eq. i35\ . which is the most 
relevant nonlinearity that violates this symmetry. This term 
can be traced back to the 72 $j term in Eq. (|2j. We define a 
dimensionless coupling constant A = X/D, where D is the 
subdiffusion constant. If we interpret the equation of motion 
as arising from a Taylor expansion of the resistivity, then we 
can write 



A 



dpo/dn 



(55) 



One expects A to be small in a not-too-dirty electron gas, since 
in a pure system pu is already non-vanishing at any nonzero 
wavevector (contributing to the denominator), while the nu- 
merator vanishes by Galilean invariance in this limit. How- 
ever, this term is strongly relevant in two dimensions. Hence, 
even if one starts with A ^ 1, in an infinite system this cou- 
pling constant will eventually become much greater than unity 
under renormalization. Ignoring this term will certainly not be 
valid in this case, so one would need to appeal to the full equa- 
tion of motion to describe the transition. In a finite system, 
however, one is interested in reducing the cutoff to roughly 
1/L, where L is the system size, so the growth of the coupling 
constant will be bounded. The model with $ — * $ + const 
symmetry will provide a reasonable description of the transi- 
tion as long as L is sufficiently small that A does not become 
of order 1. Under a tree-level renormalization group itera- 
tion, the renormalized coupling constant A' grows according 
to A' = Xs^, with s > 1. In terms of the reduced cutoff 
A' = A/s, where A is the initial cutoff, this can be expressed 
as 



X 



(56) 



We take A' = 1/Lc2 and A = 1/^in, where is the inelas- 
tic mean free path. For the samples used in the experiments, 
lin ~ hvpEp / {ksTY ~ 100/im.'^, and is comparable to the 
transport mean free path estimated from the mobility at a tem- 
perature of IK. This is about an order of magnitude smaller 
than the sample lengths. 

To estimate the critical length scale Lc2 below which the 
model is valid, we set A' = 1, leading to 



(57) 



|A| 



We note that if |A| -C 1, say around 0.01, then the criti- 
cal length Lc2 would already be comparable to the sample 
sizes studied in the experiments. A serious estimate of this 
length would require a microscopic calculation of the resis- 
tivity /0_D(k, uj) in the presence of microwaves to compute the 
bare value of A via Eq. i55\ and would be valuable. 



We have seen in the case of short-range interactions above 
that a finite-volume region of initial couplings are marginally 
irrelevant and flow to zero upon coarse-graining. Next, we 
discuss the fate of these flows when long-range interactions 
are turned on. This case is relevant experimentally due to the 
absence of metallic gates in the experiments conducted so far, 
leading to unscreened Coulomb interactions. 

Consider again the equation of motion given in Eq. i35\ . 
with ^[n] — J^, V{r — r')n(r'). With long-range Coulomb 
interactions, the Fourier-transformed interaction potential (in 
two dimensions) is V{k) ^ 1/k. The upper critical dimen- 
sion in this case is due — 3. Since we are interested in the 
transition in d = 2 dimensions, one option is to carry out an 
e-expansion in d — 3 ~ e dimensions. This approach is com- 
plicated by the need to generalize the interactions in Eq. ( I35> 
to higher dimensions. Alternatively, one can perform an e- 
expansion by writing V{k) — l/fc*^, with e ^ 1. The upper 
critical dimension is then due = 2 + e. We will adopt the latter 
approach since we can then work directly in d = 2 dimensions 
and thereby avoid generalizing the equation of motion. 

We use the dynamical renormalization group as outlined 
above to calculate the flow equations at one-loop order and 
to lowest order in e. As in the short-range case, there are no 
diagrammatic corrections to the noise strength g' at one-loop. 
To keep g' fixed under renormalization, we take the rescaling 
exponent x = (4 + z)/2. Similarly, we choose the exponent 
X ~ 3z/2 to fix the coefficient of dtn to be unity in Eq. ( I35> . 
To simplify the flow equations for the remaining parameters, 
we again use the dimensionless coupling constants defined in 
Eq. (04} (with a = g/AnD^A''). The subdiffusion constant 
then flows according to 



diD = (z-4 + e + Ai/2)D. 



(58) 



For convenience we choose z = 4 — e — l/2Ai to keep D 
fixed. 

With this choice of rescaling exponents, the flow equations 
for the parameters A^ are 



diAi 


- eAi 




+ 13A4)Ai - 3A3A4 


f 4A2A5 (59) 


diA2 


= eAa 




+ 7A4)A2 + i(3A3- 


7Ai)A5 (60) 


deA3 


= eAs 


- (-Ai 


+ 12A4)A3 


(61) 


diAi 


= eA4 


-(Ai + 


9A4)A4 + A^ 


(62) 


diA, 


= eAs 


-(Ai + 


10A4)A5. 


(63) 



In the case of short-range interactions we found that there are 
stable trajectories where all the coupling constants go asymp- 
totically to zero. This clearly cannot happen in the case of 
finite-range interactions due to the eA^ terms above. Instead, 
we search for fixed points of the form A^ = a^e, where 
are constants. One can easily show that all such fixed points 
are unstable. We interpret this lack of a stable fixed point 
as signaling a first-order transition. Thus, we conclude that 



the continuous transition that can occur with short-range in- 
teractions is driven first-order by the presence of long-range 
interactions of the form V{k) = l/fc^ with e ^ 1. 

This result may seem surprising initially since one might 
expect long-range interactions to suppress density fluctuations 
and thereby further stabilize the Gaussian fixed point. For in- 
stance, in the linearized equation of motion the density dif- 
fuses faster with long-range interactions. A competing effect, 
however, is that density fluctuations can interact nonlocally 
through the nonlinear terms. Thus, density fluctuations in one 
region of the sample can further induce fluctuations over long 
distances. This can lead to positive feedback of these density 
fluctuations via the nonlinearities, which evidently drives the 
transition first-order 
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IV. TRANSITION TO DENSITY-ORDERED STATE AT 
ZERO WAVEVECTOR IN THE FULL ROTATIONALLY 
INVARIANT MODEL 



FIG. 2: "Free energy" density / as a function of constant density n 
for A = 1, ft = 1/5, with f = 2fc, Vc, (see text for details). 



The model considered above with $ — ^ $ + const symme- 
try is only appropriate for describing physics up to a certain 
length scale. For instance, in the ordered state with a given 
uniform current, regions of negative density appear if the sam- 
ple is too large. We estimate this length scale to be roughly 
^ 1mm using parameters from Willett's experiments.^ On 
larger scales, terms in the equation of motion depending on 
the magnitude of which prevent the density from becom- 
ing negative, must be taken into account. As discussed above, 
at the transition, the leading term involving the magnitude of 
<i> (i.e., the 72<i>j term in Eq. is strongly relevant in two di- 
mensions. The dimensionless coupling constant for this term 
therefore grows under renormalization. In a finite system, 
the growth of this coupling is limited by the system size L 
since one only reduces the wavevector cutoff to of order 1/L. 
Neglecting terms depending on the magnitude of $ becomes 
an invalid approximation when the system size is sufficiently 
large that this renormalized dimensionless coupling becomes 
of order unity. 

To describe physics in samples with linear dimensions 
larger than these length scales, one must therefore relax the 
$ ^ $ + const symmetry and appeal to the full equation of 
motion in Eq. (|2} with no additional symmetries. This is the 
subject of the present section. Identifying the ordered state 
that develops in this case is nontrivial, so we will focus only 
on the transition to the ordered state, considering both short- 
and long-range interactions. 



since the diffusive mode becomes unstable when f < 0. Other 
nonlinearities are in principle present in Eq. (I64> . but they are 
less relevant than the A term and can be neglected provided 
we work near the upper critical dimension. 

To derive Eq. ( I64> . we solved for the current in terms of 
the density assuming two spatial dimensions. However, the A 
interaction is strongly relevant in two dimensions, so to study 
its effects we need to to continue this model to higher dimen- 
sions. We initially adopt the most naive way of doing this con- 
tinuation; namely, we simply assert that Eq. (I64> is valid in d 
dimensions. In the case of short-range interactions the upper- 
critical dimension for the A nonlinearity is then d^c = 6. 

We have carried out an e-expansion in d = 6 — e dimen- 
sions to obtain the renormalization group flow equations to 
one-loop order Rather than go through the details of the cal- 
culation, we will merely state that these equations lack a stable 
fixed point, which we interpret as signaling a first-order transi- 
tion. A more direct route to this conclusion can be obtained by 
observing that Eq. ( I64> is identical to the equation of motion 
for an equilibrium model with a conserved order parameter. 
That is, it (Eq. (I64» can be rewritten as 



d,n = V^'-f + V-C', 
on 



with the "free energy" given by 



F ■ 



fn^ + D{Vnf 



(65) 



(66) 



The term proportional to u that results from Eq. ( I65> is irrel- 
evant at the upper critical dimension and has therefore been 
excluded from Eq. (I64> . We will assume u > for simplicity, 
although this is not essential. Figure|2depicts the free energy 
density / as a function of uniform density n for three differ- 
= dtn - fV^$ + DV^^ - V • C' - AV ■ ((1>V$), (64) values of f. When f > A^/18u ee f^, the free energy is 

minimized when n = as illustrated by the solid curve. At 
f — fc, the free energy has two degenerate minima as shown 
in the dashed curve. Below fc, the free energy is minimized 
by a nonzero value of n. This situation is represented by the 



A. Transition with sliort-range interactions 

When we relax the $ ^ $ + const symmetry, Eq. ( I35> 
generalizes to 



where D > and is a Gaussian noise source with variance 
2g'. In this subsection we consider short-range interactions, so 
= n. The transition in the linearized theory occurs at 7~ = 
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dotted line for f = 0. When f decreases below rc, the density 
will therefore jump discontinuously from zero to minimize the 
"free energy". This signals the onset of a first-order transi- 
tion, consistent with our renormalization group results. Note 
that the transition occurs at a finite value of r, preempting the 
apparent transition (a 'spinodal') at f = expected from the 
linear theory. As Fig. |2] demonstrates, the point f = actu- 
ally corresponds to a spinodal decomposition where the sys- 
tem goes from being metastable to globally unstable at n = 0. 

These results hold only near d = 6. We can reduce the up- 
per critical dimension of the model by considering a spatially 
anisotropic continuation of Eq. i64\ to higher dimensions. To 
do this, we can start by continuing Eq. (|2j to d dimensions and 
taking the resistance at zero wavevector and frequency to be 
anisotropic. That is, write 

rj ^ r±j±+ri\j\\ (67) 

in Eq. (|2j, where j± is the current in the x-y plane and j || rep- 
resents the current in the additional d — 2 dimensions. We will 
be interested in tuning the resistance r± in the x-y plane to 
zero while leaving the resistance r|| for the remaining direc- 
tions positive. We can then eliminate the current in favor of 
the density as before to obtain 

= atn-f||V|$ + i:'±vl$ 

- Vi-Cl-AiVi-($Vx$), (68) 

where f\\ , D± > and we have set rj_ = 0. We have also 
only retained the x-y (in-plane) components of the noise <^'j_, 
as noise components in the additional d — 2 || dimensions 
are irrelevant. Similarly, we have omitted nonlinear terms in- 
volving V|| since, due to the high anisotropy of the harmonic 
terms, these are clearly less relevant than the corresponding 
terms involving only \7± derivatives. 

The upper critical dimension for this model is due = 4. We 
have performed an e-expansion in d = 4 — e dimensions to 
one-loop order Once again, we find that the model lacks a 
stable fixed point. Thus, even near d = 4 dimensions, the 
transition still appears to be first-order 

It seems quite likely that the transition is first-order in d = 2 
dimensions as well. In the model with $ ^ $ + const sym- 
metry and short-range interactions, we showed in the previous 
section that one could have a continuous transition in d = 2 
dimensions. Terms that violate the $ ^ $ + const symmetry 
appear then to always drive the transition first-order. 

We propose the following simple physical interpretation for 
this. We have been analyzing the transition at zero wavevector 
and zero frequency, where one expects long-wavelength fluc- 
tuations that become critical at the transition to give rise to 
an ordered state with uniform, static current. Such an ordered 
state must be accompanied by a density gradient transverse to 
the current to balance the Lorentz force. We have already ar- 
gued that such a state cannot exist in the thermodynamic limit 
because the density would become arbitrarily large and nega- 
tive at the edges of the sample. The only terms in the equation 
of motion that sense these unphysical features are precisely 
those terms that depend on the magnitude of the density. In 
the thermodynamic limit, these terms must therefore induce 



a first-order transition into some other state, such as a state 
ordered at finite wavevector or a phase-separated state. A di- 
rect transition from a uniform isotropic liquid to a modulated 
(finite wavevector) smectic state can also be argued to be first- 
order on quite general groundsi2i2^i2^ 

B. Transition with long-range interactions 

Finally, let us consider the effect of long-range interactions. 
We saw in the model with $ $ + const symmetry that 
turning on long-range interactions drove the transition first- 
order. In the present case, the transition is already first-order 
with short-range interactions, so it seems rather likely that the 
transition will remain so with long-range interactions. This is 
indeed what we find based on a renormalization group anal- 
ysis. We will therefore only outline the calculation and state 
the results. 

Consider Eq. i64l with ^[n] — J^, V{r — r')n(r') and 
V{k) = 1/fc. Once again, the A interaction is strongly rele- 
vant in two dimensions, so we would like to continue Eq. i64l 
to d dimensions. We will only consider the simplest isotropic 
continuation and assert that Eq. ( I64l i holds in d dimensions. 
The upper critical dimension is then d^c — 7. We have per- 
formed a one-loop e-expansion in d = 7 — e dimensions, and 
find that the model lacks a stable fixed point. Thus, as ex- 
pected, the transition remains first-order when long-range in- 
teractions are included. 



V. DISCUSSION AND SUMMARY 

The focus of this paper has been on the physics near the 
transition to zero resistance state in 2DEGs driven with mi- 
crowave radiation. Our goal was to understand the long- 
distance, long-time properties of the system taking into ac- 
count noise and fluctuation effects within a nonequilibrium 
hydrodynamic theory involving the electron current and den- 
sity. We specifically focused on the transition to a time- 
independent, density-ordered state that occurs when the mi- 
croscopic resistance first becomes negative at (k = 0, w = 0). 
The long wavelength subdiffusive density fluctuations are the 
only critical modes at this transition. We analyzed two models 
involving the density mode: (i) Model-1, characterized by an 
imposed symmetry under a global uniform shift of the density, 
valid only on sufficiently small length scales, and (ii) Model- 
11, which is most general rotationally invariant model with no 
additional symmetries. 

The ordered state in Model-1 consists of a uniform current 
and a transverse Hall electric field that balances the Lorentz 
force. This state was shown to be stable within a linearized 
theory of fluctuations about the ordered state. 

We argued that the uniform-current steady state in Model-I 
cannot exist in arbitrary large samples since the uniform Hall 
field would eventually lead to regions of negative density. Us- 
ing parameters from Willett's experiments^, we estimated that 
samples with dimension smaller Ld ^ 1mm can support this 
state. To describe the ordered state in larger samples, one must 
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include terms that depend on the magnitude of the density, 
which would prevent the density from becoming arbitrarily 
large and negative. 

Since the ordered state in Model-1 breaks continuous rota- 
tional symmetry, there is an associated Goldstone mode corre- 
sponding to long-wavelength fluctuations of the current trans- 
verse to the uniform current flow direction. We suggest that 
a possible way of detecting this Goldstone mode might be to 
use surface acoustic waves in the zero resistance regimeiS2& 
A surface acoustic wave at the right wavelength and frequency 
should couple to this excitation, leading to anomalous shifts in 
the velocity and intensity of the wave. 

The transition to this ordered state in Model-I was analyzed 
both in the case of short- and long-range interactions using dy- 
namical renormalization group methods. This model is valid 
for describing the transition on length scales L < Lc2, which 
we think could be comparable to sample sizes in current ex- 
periments as discussed in Section 111 C, although it would be 
valuable to have an estimate from microscopic calculations. 
With long-range interactions, we showed that Model-1 under- 
goes a first-order transition. However, with short-range in- 
teractions, we showed that in two dimensions the Gaussian 
fixed point in Model-I has a finite-volume basin of attraction. 
That is, a finite-volume region of initial nonlinear couplings 
all flow to zero upon renormalization. The transition in these 
cases is of the mean field type. In particular, mean-field the- 
ory predicts a continuous transition to the ordered state. Ad- 
ditionally, the density subdiffuses at the critical point, with a 
frequency given by w cx — i/c"*. This subdiffusion should lead 
to large density fluctuations and hence large voltage fluctua- 
tions at the transition. It may be interesting to observe this in 
samples with metallic gates, so that the Coulomb interactions 
are screened. Although it may be difficult to quantitatively 
test the mean-field predictions, one could perhaps measure 
voltage correlations at contacts placed along the perimeter of 
the sample. These voltages should behave similarly to the 
density-density correlations given in Eq. i52\ . Qualitatively, 
one should at least observe large voltage fluctuations since the 
density is critical and subdiffusive at the transition. 

We next turned to an analysis of the transition in the more 
generic Model-11, which includes terms that depend explicitly 
on the magnitude of the density. We found that the transition 
within this model is always first-order independent of whether 
interactions are short- or long-range, at least near the upper- 
critical dimension of the theory. The physical mechanism for 
this first-order transition is as follows. If the resistance min- 
imum occurs at fc = w = 0, then one would expect long- 
wavelength fluctuations to give to a time-independent state 
with a uniform current and transverse Hall field. As men- 
tioned above, such a state cannot exist in arbitrarily large sam- 
ples since regions of negative density would eventually ap- 
pear. The role of terms that depend on the magnitude of the 
density is to prevent such unphysical features from arising. 
These terms consequently force a first-order transition into a 
more complicated ordered state. 

The experiments conducted so far were carried out us- 
ing samples without metallic gates, leading to unscreened 
Coulomb interactions. The transition in these systems is there- 



fore predicted to be first-order, which should have measur- 
able consequences. In particular, one would expect discon- 
tinuous jumps in various observable quantities such as spon- 
taneous currents, voltages, and local magnetizations that de- 
velop at the transition. We realize that these jumps may be 
difficult to measure experimentally. Another possibly more 
controlled way of detecting a first-order signature might be to 
measure the critical current above which the zero-resistance 
state disappears.— If one approaches the transition from the 
ordered state (by, say, changing the magnetic field) then the 
critical current should drop discontinuously from some finite 
value to zero if the transition is indeed first-order. 

There are several future directions one could pursue with 
the theory presented here that we believe would be interest- 
ing and provide further insight into the remarkable physics of 
driven 2DEGs. Regarding the transition to zero resistance, we 
have considered only the simplest case where the resistance 
minimum occurs at fc = = 0. It may be interesting to gen- 
eralize our results for this case to include static disorder to see 
how it affects the transition. One could also analyze the tran- 
sition at finite frequency where a time-dependent state such as 
circulating currents would arise. Additionally, one could con- 
sider the transition at zero frequency but nonzero wavevector 
fco. In this case one would be interested in wavevectors k such 
that |fc — fcol < A for some cutoff A. If the equation of motion 
was derivable from a free energy of the form 

F ^ r(q)n(q)n(-q) 

+ / A(qi,q2,q3)n(qi)n(q2)n(q3) 

+ / Ii(qi,...,q4)n(qi)n(q2)n(q3)»^(q4)(69) 

J qi- ■q4 

then we know that the cubic term drives the transition first- 
order based on analogies with the solidification of an isotropic 
liquidi2i Even in the case where the cubic term vanishes, the 
transition is still driven first-order by fluctuations^^. Due to 
the presence of nonequilibrium terms, however, the equation 
of motion will not be derivable from a free energy. Nonequi- 
librium effects could cause dramatic deviations from the equi- 
librium theory, and at present it is unclear what effect such 
terms will have on the transition. 

Another avenue one could pursue with this theory is to ad- 
dress the properties of the ordered state away from the transi- 
tion in Model-II. Numerical studies may be best suited for this 
purpose especially since the ordered state is likely to be inho- 
mogeneous and not analytically tractable. One possible route 
is to generalize the numerics on the flocking transition done 
by Vicsek et al}'^ to include a magnetic field and interactions. 

Finally, we note that while early numerical work on the 
flocking transition in the absence of a magnetic field indicated 
a continuous phase transition^^, some recent simulations on 
larger system sizes hint at a weak first-order transition^". If 
true, this would be consistent with the transition continuing 
to be first-order in the presence of a magnetic field for large 
enough systems as argued in this paper. 
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